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ABSTRACT 

We performed micro-gravity collision experiments in our laboratory drop-tower using 5-cm-sized 
dust agglomerates with volume filling fac tors of 0.3 a n d 0.4 , respectively. This work is an extension 
of our previous experiments reported in iBeitz et al.l (|2011l) to aggregates of more than one order 
of magnitude higher masses. The dust aggregates consisted of micrometer-sized silica particles and 
were macroscopically homogeneous. We measured the coefficient of restitution for collision velocities 
ranging from 1 cm s _1 to 0.5 m s -1 , and determined the fragmentation velocity. For low ve locities, the 
coeffic ient of restitution decreases with increasing impact velocity, in contrast to findings bv Be itz et al.l 
(|2011f ). At higher velocities, the value of the coefficient of restitution becomes constant, before the 
aggregates break at the onset of fragmentation. We interpret the qualitative change in the coefficient of 
restitution as the transition from a solid-body-dominated to a granular-medium-dominated behavior. 
We complement our experiments by molecular dynamics simulations of porous aggregates and obtain 
a reasonable match to the experimental data. We discuss the importance of our experiments for 
protoplanetary disks, debris disks, and planetary rings. The work is an extensional study to previous 
work of our group and gives a new insight in the velocity dependency of the coefficient of restitution 
due to improved measurements, better statistics and a theoretical approach. 

Subject headings: coefficient of restitution — fragmentation — dust — planetary rings — debris disks 
— impact — solar system: formation 



1. INTRODUCTION 

There is observational evidence that cm-sized par- 
ticle s exist in protop lanetary disks (PPDs) (see 
e.g. iWilner et al.l 120051) . A recent model of the 
protoplanetary-dust growth, based on laboratory ex- 
periments, has shown that cm-siz ed particles can be 
formed by direct c ollisional sticking (IGiittler et al J 1201 01 : 
IZsom et all 120101 ) . IGiittler et all (|20100 also found that 
collisional compaction can lead to filling factors of up to 
0.57. Compaction in fra gmenting coll i sions with mass 
transfer was also fou nd bvlKothe et al.l (120 101. who con- 
firmed the model bv IGiittler et al.l (|2010|^ 

If there were even larger solid particles avail- 
able in PPDs, then the growth could commence 
through a fragmentation-coagulation cycle, leading to 
dust-aggregates sizes in the planetesimal size range 
(|Windmark et al.ll2012t ). As the dust growth in PPDs 
starts with (sub-)/im dust grains, it is natural to assume 
that the resulting macroscopic bodies are agglomerates 
of the microscopically small dust grains. Many bodies 
in debris disks and planetary rings are also expected to 
consist of such granular material. Therefore, it is inter- 
esting and important to know the collision behavior of 
very large dust aggregates. A first approach to the >cm 
size range was established bv IBeitz et al.1 (j2011[ ). who in- 
vestigated collisions among 2-cm-sized dust aggregates. 
Here, we present follow-up experiments with 5-cm-sized 
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dust aggregates, which are more than an order of mag- 
nitude larger in mass. 

In Section [2] and E2 we describe our experimental ap- 
proach and the experimental results. Section |4] explains 
the numerical model to understand the physics in dust- 
aggregate collisions, and Section [5] gives some astrophys- 
ical applications for low-velocity collisions of granular 
bodies. Finally, Section [6] summarizes our results. 

2. EXPERIMENTAL APPROACH 

In this Section, we describe the experimental methods 
applied for the determination of the coefficient of restitu- 
tion and the fragmentation threshold of 5-cm-sized dust 
agglomerates. 

2.1. Preparation of the Dust Agglomerates 

The dust material for the production of large aggre- 
gates was pure Si02 powder, consisting of 0.1-10/im sized 
irregular grains. Owing to the fact that direct growth can 
commence through the mm-size range, upon which new 
growth process (e.g., fragmentation with mass transfer) 
occur, we expect that large dust agglomerates in PPDs 
possess a hierarchic structure and are agglomerates of 
agglomerates. Therefore, prior to the compression into 
large dust aggregates, the dust powder was being sieved 
through a mesh with a width of 0.5 mm to av oid larger 
aggregates. As shown bv lWeidling et al.l (|2012l ). the siev- 
ing process produces dust agglomerates with a filling fac- 
tor of 0.35. The sieved dust was then filled into a hollow 
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steel cylinder with 5 cm diameter and was then slowly 
compressed with a brass piston. The volume filling fac- 
tor of the compressed dust aggregates was adjusted by 
filling the cylinder with a defined dust mass m and by 
compressing the sample to a pre-determined height of 5 
cm. After the compression, the bottom of the steel cylin- 
der was removed and the dust agglomerate was pushed 
out of the cylinder using the piston. 

2.2. Morphology of the Dust Agglomerates 

To examine the inner structure of our dust agglom- 
erates, we performed x-ray tomography (XRT) measure- 
ments. Figure [1] shows a plane perpendicular to the cylin- 
der axis of 5 cm diameter and shows granularity, which 
represents the original sieved dust aggregates. This sub- 
structure is not visible in some areas of a few mm in 
radius. We think that this is caused by slightly inho- 
mogeneous compression. To increase the resolution of 
the XRT, we cut a mm-sized piece out of a larger dust 
aggregate, which included the cylinder mantle, and re- 
peated the XRT measurements with a voxel size of 17/im 
x 17^m x 17/um. In the inset of Figure [TJ a cut through 
the center of the sample is displayed, with the former 
cylinder mantle at the top of the inset picture. On the 
cylinder mantle, the agglomerate-of-agglomerates struc- 
ture is not present. The material has obviously formed 
a slightly densified mantle of about 85/im thickness with 
a filling factor of 0.35. We do not think that this slight 
density increase by about 10 percent has a considerable 
effect on the impact behavior of the dust agglomerates. 




Fig. 1. — Reconstruction of an XRT image showing a radial slice 
perpendicular to the cylinder axis in the center of the dust agglom- 
erate (window size: 5.9 cm X 5.9 cm). Inset: cut-out of a piece of 
the dust agglomerate, including the cylinder mantle (window size: 
0.63 cm X 0.54 cm). 

Figure [2] shows the density distribution inside the dust 
agglomerate. In the top graph, the mean density of 
planes perpendicular to the cylinder axis is shown over 
the full height of the agglomerate. The median volume 
filling factor of the dust agglomerate is ~ 0.3. At the top 
6 mm of the dust-agglomerate cylinder, the mean density 
is increased to a filling factor of up to 0.38. This is the 



place where the piston has pushed against the agglomer- 
ate during compression. In the lower 5mm of the dust ag- 
gregate, the density of the agglomerate rapidly decreases 
to a volume filling factor of about 0.15. The lower graph 
of Figure [5] shows the radial density distribution averaged 
over the azimuth and over 20% of the cylinder height for 
each curve. The upper curve was averaged over the up- 
per 20% of the dust agglomerate. The three curves in 
the center were averaged over the central 20%-40%, 40- 
60% and 60%-80% slices of the dust agglomerate. Mind 
that these three curves have basically the same values 
and overlap in Figure [5J The lower curve was averaged 
over the lower 20% of the dust-agglomerates height. The 
lower graph of Figure [2] shows that the density of the 
central part of the dust agglomerate is slightly increasing 
from the axis to the mantle by about 5% in volume filling 
factor. Here, the top and bottom curves are of no partic- 
ular importance, because the dust agglomerates collide 
close to their midplanes. Within the accuracy of the 
XRT measurements (~ 10— 100/im), only slightly densi- 
fied material on the dust-agglomerate's cylinder mantle is 
found, which we believe does not dramatically influence 
the collisional outcome. 

2.3. Collision Experiments 

The collision experiments using the dust aggregates de- 
scribed above are performed in our 1. 5-m laboratory drop 
tower described in lBeitz et al.l (|201lD . The agglomerates 
were placed above each other with their symmetry axes 
rotated by 90 degrees (see Figure [3] top) . Each dust ag- 
glomerate is supported by two brackets, which can be 
rapidly pulled away by solenoid magnets. The upper dust 
agglomerate is released slightly earlier than the lower one 
(see Figure [3]). The time difference At between the re- 
lease of the upper and lower dust agglomerate results in a 
relative velocity v — gAt between the two bodies, where 
g = 9.8 m s~ 2 is the surface acceleration of the Earth. 
As the drop tower is evacuated to a residual gas pressure 
of 100 Pa, air drag can be neglected during the ~ 0.5 s 
free-fall time. 

The resulting impacts are observed by two cameras 
(one high-speed camera, one camera with a normal frame 
rate) with 90 degrees angular distance. The cameras, 
which are outside the drop tower, are released such that 
they fall in the center-of-mass frame of the two dust ag- 
gregates. 

The release mechanisms of the two dust agglomerates 
cause a slight rotation of the two bodies. Furthermore, 
the collisions were not always perfectly central, which 
also causes rotation. To account for these effects, we 
calculated the coefficient of restitution from the ratio of 
the square root of the total kinetic energies of the two 
dust agglomerates in their center-of-mass frame before 
and after the collision, 

= J m *i + ff + "f \. (i) 

Here, V\ and i>2 are the relative translational velocities of 
the two dust aggregates before (index 1) and after (index 
2) the collision, w Ul W12, W21, and W22 are the mutually 
perpendicular angular velocity components of the collid- 
ing cylinders (the first index refers to before and after 
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Fig. 2. — Volume filling factors derived from the XRT measure- 
ments of a dust agglomerate. Top graph: the mean density of 
planes perpendicular to the cylinder axis as a function of the height 
of the dust agglomerate. Bottom graph: the radial density distri- 
bution averaged over the azimuth and 20% of the dust-cylinder 
height for five vertical positions of the dust agglomerate. The up- 
permost curve belongs to the top part of the dust aggregate, the 
three center curves (partly overlapping) belong to its central parts 
and the lower curve represents the lower portion of the dust ag- 
glomerate. 



the collision, the second index describes, which compo- 
nent of the angular velocity is meant), and m* = m/2 
and I are the reduced mass and the moment of inertia 
of the two dust aggregates. Please note, that we only 
consider rotation of the dust cylinders around axes per- 
pendicular to their symmetry axis. A rolling motion of 
the dust aggregates, which is both difficult to observe and 
possesses a much smaller moment of inertia, is neglected 
here. The effective relative velocity, including rotation 
effects, of the agglomerates prior to the collision is then 
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Fig. 3. — Top: the two dust agglomerates shortly after release 
and before the collision. Both release mechanism are fully open. 
Bottom: the two dust aggregates fragment shortly after impact. 

given by 

l m*vf + I(Lof 1+ ^ 12 ) 

y ef£ = V ; • ( 2 ) 

V m* 

This rotation of the dust aggregates is only important 
at the very lowest impact velocities; at higher relative 
velocities, the rotational motion prior to the impact is 
negligible with respect to the rotation after the impact. 

At velocities lower than <~ 10 cm s^ 1 , the statistical 
scattering of the single measurement was very high. This 
was due to the fact that the free-fall time of the dust 
aggregates was limited to ^0.5 s, which allows for a 1 
cm s^ 1 impact velocity a maximum distance between the 
dust agglomerates prior and after the impact of only 2.5 
mm. This is both from the preparatory point of view and 
for data-analysis reasons a limiting value so that lower 
impact velocities were not achievable. Therefore, more 
experiments have been performed at the lower velocities 
to get a reasonable mean value for the coefficient of resti- 
tution. 

3. EXPERIMENTAL RESULTS 
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In this Section, we will present our data on the low- 
velocity coefficients of restitution as well as on the frag- 
mentation velocities of dust aggregates of two different 
volume filling factors. 

3.1. Coefficient of Restitution 

In Figure 21 the coefficient of restitution of our 5-cm 
sized cylindrical agglomerates, according to the definition 
in Eq. [1] is plotted as a function of the effective impact 
velocity, according to Eq. [2] for volume filling factors of 
4> = 0.3 and <fi — 0.4, respectively. Each data point is the 
mean of four measurements; the error bars denote the 
2-cr error of the mean value. 
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Fig. 4. — The coefficient of restitution as a function of the ef- 
fective impact velocity for 5-cm sized dust agglomerates of 30% 
volume filling factor (upper graph) and 40% volume filling factor 
(lower graph), respectively. Each data point is the mean value of 
four measurements; the error bars denote the 2-a uncertainties of 
the mean values. The dotted line shows a power law with a s l ope of 
-1/4, follow ing t he solid-state model by [Thornton & Nina H1998I ) 
(sec Section l4.ll) . 



We can clearly distinguish different velocity regions: 
(1) At the very lowest velocities (v e g < 4 cm s _1 ), the 
coefficient of resti t ution follows the solid-state model by 
iThornton fc Ningl (| 19981 ) within the measurement uncer- 
tainties as shown by the error bars in Figure @] (see Sec- 
tion [4~T]). (2) At higher velocities (4 cm s _1 < v e g < 



6 cm s _1 ) for a filling factor of 0.3 and 4 cm s" 1 < 
v ff 10 cm s 1 f° r a filling factor of 0.4), the coefficient 
of restitution decreases steeply to values of e — 0.2 for 
cj) = 0.3 and e = 0.3 for <j> — 0.4, respectively. (4) At even 
higher velocities, we observe fragmentation (see Section 
1331) . 

The coefficient of restitution was measured by consid- 
ering translational and rotational velocity changes. Be- 
cause translational and rotational velocities are trans- 
formed into one another, depending on the impact pa- 
rameter, it is not possible to split the coefficient of resti- 
tution into components of rotation and translation. To 
show the contribution of rotation in our measurements, 
we display in Figure [5] the ratio of the translational ve- 
locity to the effective velocity prior to the impact as well 
as the ratio of the coefficients of restitution with and 
without the consideration of rotation over the effective 
impac t vel ocity We use the same averaging process as 
in fig. 14.11 Figure [5] shows that the rotational contribu- 
tion in our measurements is below 40% over all veloci- 
ties and that the influence of rotation on the coefficient 
of restitution increases with decreasing impact velocity. 
However, the data also show that the different regimes of 
the coefficient of restitution (see above) are not caused 
by rotation effects. 

In contrast to our r esult, iBeitz et al.l (|2011h found 
no obvious correlation between the coefficient of resti- 
tution and impact velocity. Unlike in our experiments, 
they used smaller spherical dust agglomerates made from 
monodisperse spherical grains. We used cylindrical ag- 
glomerates, because we were not able to produce homoge- 
neous crack-free spherical dust agglomerates. Our poly- 
disperse irregular SiC>2 grains possess smaller contact 
forces between the mon omer particles tha n the monodis- 
perse spherical grains of Be itz et all ([201 If) . This and the 
smaller filling factor in our experiments ease collisional 
compression and reduce the effect of randomly jammed 
packi ngs during collision (see e.g. iTorquato fc Stillingerl 
I2010H . Jammed packings are stronger at higher compres- 
sion velocities and increase the rigidity of an agglomerate 
and therewith decrease the contact ar ea of the colliding 
agglom erates. Due to the model by IThornton fc Nina 
(1998), this should increase the coefficient of restitution. 
Therefore, it is possible th a t jam ming occurred in the 
experiments by IBeitz et al.l (|2011| ) , due to their higher 
filling factor and monodisperse spherical grains, which is 
the more pronounced at higher the collisional velocities. 

3.2. Fragmentation Velocity 

Above a certain impact velocity, one or both of the 
colliding dust aggregates fragment so that the determi- 
nation of the coefficient of restitution becomes meaning- 
less. We measured the outcome of collisions among the 
5-cm sized dust aggregates with volume filling factors 
of <j> = 0.3, <j) — 0-4 and </> = 0.48 for velocities up to 
200 cm s _1 and plotted the collision results in Figure [51 
The + signs denote bouncing, A stands for the fragmen- 
tation of one of the two dust aggregates, whereas * and o 
describe the fragmentation of both dusty bodies (in the 
latter case with a maximum fragment mass smaller than 
half of the pr ojectile mass ) . For comparison, we added 
the data from feeitz et alj ([201 If ) for 2-cm sized spheri- 
cal dust aggregates with 4> — 0.5, using spherical 1.5^m 
sized SiC>2 monomers. Within the scatter of data, the 
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Fig. 5. — The ratio of the effective reduced velocity (including 
rotation) to linear velocity prior to impact (+) and the ratio or the 
coefficient of restitution including rotational effects to the coeffi- 
cient of restitution with neglected rotation (*). 

threshold between bouncing and single- aggregate frag- 
mentation seems to be independent of the volume filling 
factor at Uf rag = 40 ± 10 cm s _1 . This is about a fac- 
tor oftwjoJngher_Hmn the fragmentation threshold found 
bv lBeitz et al.l (|2011[) and might be due to the different 
geometry (spheres vs. cylinders), the different monomer 
morphologies and sizes (irregular polydisperse vs. spher- 
ical monodisperse), dust-aggregate size (5 cm vs. 2 cm) 
or volume filling factor (<fi < 0.48 vs. 4> — 0.5). 

4. MODELING DUST-AGGREGATE COLLISIONS 

4.1. Continuum Theories 

Assuming that at low velocities dust agglomerates 
behave like solid bodie s, we can apply the theory of 
iThornton fc Ningl (|1998|) for the derivation of the impact- 
velocity dependence of the coefficient of restitution. At 
larger velocities, we assume that the granular material 
gets more and more mobilized and begins to fluidize so 
that a so lid-state theory is no lo nger applicable. The 
model by IThornton fc Ninel (|1998f) relies on two thresh- 
old velocities, the transition velocity between sticking 
and bouncing and the transition velocity between elas- 
tic and plastic material effects, v y , respectively. The 
data shown in Figure [4] suggest that our collision ve- 
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Fig. 6. — The possible outcomes in collisions among 5-cm sized 
dust aggregates for different volume filling factors between <f> = 0.3 
and cf> = 0.48. A + denotes bouncing, A stands for the fragmen- 
tation of one of the two dust aggregates, * and o describe the 
fragmentation of both dusty bodies, the latter with a maximum 
fragment mass smaller than half o f the p ro jectile mass. For com- 
parison, the data from Bcitz ct al. (2011) for 2-cm sized spherical 
dust aggregates with <j> = 0.5, using spherical 1.5/^m sized Si02 
monomers, are shown additionally. 

locities are clearly above both thresholds, because the 
coefficient of restitution is in the regime in which energy- 
dissipating effects domi nate. In this velocity regime, 
IThornton fc Nine) (|1998l ) predict a decrease of the coeffi- 
cient of restitution with increasing velocity and asymp- 
totically a power-law behavior of the coefficient of resti- 
tution of the form e oc i> -1 / 4 . The dotted line in Fig- 
ure |4] shows such a dependence, with v y — 0.12 cm s _1 
for (j> — 0.3 and v y = 0.4 cm s _1 for <f> = 0.4, respec- 
tively. As can be seen in the comparison between the 
model and our data, our 5-cm sized dust aggrega t es can 
be reasonably described by the IThornton fc Ningl (1998) 
continuum theory for effective impact velocities up to 
~ 4 cm s _1 . However, for larger impact speeds, the co- 
effi cient of restitution drops below the curve predicted 
by IThornton fc Ninel fl998) an d is rather velocity in- 
dependent for velocities between ~ 10 cm s" 1 and the 
fragmentation limit (see Figure 0]). 

4.2. Molecular Dynamics Simulations of Dust-Aggregate 

Collisions 

To understand these different collision behaviors of 
dust aggregates, we use a molecular dynamics ap- 
proach featuring detailed contact mechanics of micro- 
scopic silicate grains. The corresponding intera ction 
laws have been proposed bv Uohnson et al.l (|1971h and 
IDominik fc Tielensl (119951 11996ft . An overview is given 
bv IDominik fc Tielensl (|1997t ) who applied the model to 
simulate dust agglomerate collisions for t he first time 
in the context of planet formation. Later, iWada et al.l 
(|2007f ) presented a different approach where the same in- 
teraction laws were d erived from potentials. Recently, 
ISeizinger et al.l (|2012f ) proposed simple modifications to 
the model to better reproduce the compr ession behav- 
ior me asured in laboratory experiments bv lGuttler et all 

<M)'. , 
In t his work the modified model of Seiz inger et al.l 

(2012) is used. The material parameters used here 



6 



Schrapler, Blum, Seizinger & Kley 



TABLE 1 

Material Parameters of the individual monomers used in 
the molecular dynamics simulations. 



Physical property Silicate 



Particle Radius r (in /im) 0.6 

Density p (in g cm -3 ) 2.65 

Surface Energy 7 (in mj m — 2 ) 20 

Young's Modulus E (in GPa) 54 

Poisson Number v 0.17 



Critical Rolling Length £ cr it (in nm) 2 




Fig. 7. — An example of an aggregate used for the collision sim- 
ulations. It possesses a mean filling factor of <j> = 0.55 and consists 
of pa 40, 000 monomers. 

are given in Table [T] They are identical to those of 
ISeizinger et al.l (|2012l ) but differ from the ones used for 
similar simulations performed bv lWada et al.l (|201 lh ■ 

4.2.1. Sample Generation 

Here, we study the collisions of spherical aggregates. 
The samples were generated using the following proce- 
dure: we began with the regular lattice of the hexagonal 
closest packing (hep) for which the volume filling fac- 
tor is given by <j) = tt/2-\/3 « 0.74. To achieve lower 
filling factors, we removed the proper amount of ran- 
domly selected monomers. Especially for lower filling 
factors, some monomers may end up disconnected from 
the rest of the aggregate. In the last step, we therefore 
removed all isolated monomers which were not connected 
to the main aggregate (typically less than 10 -3 of the to- 
tal number of monomers). Except for the removal of iso- 
lat ed monomer s , this preparation method has been used 
bv lWada et all (|20H before. 

The samples we used in this work consist of approxi- 
mately 4 x 10 4 monomers, with a coordination number of 

8.6, a diameter of 50 /^m and a volume filling factor of 
<j> ~ 0.55. An image of such a sample is shown in Figure 

El 

4.2.2. Numerical Collision Experiments 



Three samples generated by the procedure described 
above were rotated randomly and collided head-on with 
impact velocities between 0.01 and 0.5 ms" 1 . In case 
of a bouncing event, we determined the relative velocity 
between the collision partners by averaging over the in- 
dividual velocities of the monomers that each aggregate 
was composed of. 

Due to the lattice structure of the samples, their rel- 
ative orientation has a significant influence on the out- 
come of the collision. Depending on this orientation, 
both sticking and bouncing may occur for a given veloc- 
ity. As we intended to study the dependency of the coef- 
ficient of restitution on the collision velocity, we kept the 
same orientation for a whole collision sequence. For each 
sample, eight randomly chosen orientations were exam- 
ined. Thus, 24 collisions sequences have been simulated 
in total. Typically, not all collisions of one sequence re- 
sulted in bouncing. 

4.2.3. Results 

The results of these collision sequences for the coef- 
ficient of restitution are shown in Figure [8] Similar to 
the laboratory experiments, one can distinguish different 
velocity regimes. The velocity dependency of the coef- 
ficient of restitution can be well fitted by power laws, 
e(v) — av b . We determined two fits ei ow (?;) for the low- 
velocity regime ranging from 1cm s" 1 to 7cm s -1 and 
ehigh('f) for velocities from 5cm s~ to 33cm s -1 . We 
obtain ai ow = 0.219, &i ow = —0.268 and a n igh = 0.109, 
frhigh = —0.513. For the low-velocity regime, the ex- 
ponent of the power law agrees ver y well with the one 
derived by iThornton &; Nina ([1998) for the continuum 
theory, i.e. ^theory — —0.25. The error bars show the 
influence of the orientation of the colliding aggregates on 
the outcome of the collision. For velocities higher than 
33cm s _1 , we did not observe bouncing behavior any- 
more. 
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Fig. 8. — Dependency of the coefficient of restitution on the col- 
lision velocity obtained from molecular dynamics simulations. The 
collision partners are two identical £s 50 /an sized spherical dust 
aggregates having a mean porosity of <f> pa 0.55. 

For low-velocity collisions, the kinetic impact energy is 
too low to restructure the aggregates and, thus, the en- 
ergy is mainly dissipated by the formation and breaking 
of contacts between the monomers as shown in Figure [9j 
At velocities above 2 cm s _1 , restructuring of the dust 
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aggregate sets in. As we can see in Figure [TU1 inelastic 
rolling and sliding become the most important dissipative 
channels. The dominance of the inelastic sliding results 
from the high compactness of the aggregates. Due to the 
large coordination number, monomers are tightly fixed 
by their neighbors, which limits the amount of inelastic 
rolling that may occur. 
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Fig. 9. — Energy dissipation during a bouncing collision in the 
low velocity regime with an impact velocity of 1 cm s" 1 . The 
impact occurs after 23 fis. The kinetic energy is dissipated by 
contact formation and breaking and vibrations of monomers. 
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Fig. 10. — Energy dissipation of a bouncing collision in the high 
velocity regime with an impact velocity of 10 cm s . The impact 
occurs after 2.3 fim. Contrary to the low-velocity regime (Figure 
(9{, energy is mainly being dissipated by inelastic sliding and rolling. 

Given the vast differences in monomer-particle num- 
bers between the experimental and the model dust ag- 
gregates, the agreement in the range of absolute values of 
the coefficient of restitution and the occurrence of differ- 
ent physical regimes is striking. Both approaches show 
a solid-state-like behavior for very low impact velocities 
and deviations from the expected power law with a slope 
of -1/4 for higher velocities. 

5. ASTROPHYSICAL APPLICATIONS 
5.1. Protoplanetary Disks 



It has only recently been sh own that cm-size d dust 
aggregates can grow in PPDs ()Zsom etaHl2010h . Un- 
der certain co nditions, even lar g er du st aggregates can 
be formed. iWindmark et al.l (|2012| ) showed that a 
few indestructible, cm sized solid bodies can trigger 
the further growth of dust aggregates through a frag- 
mentation and re-accretion cycle. Baroclinic vortices 
( Klahr fc Bodenheimerl I2004D or streaming instabilities 
(lYoudin fc Goodman! 120051 can concentrate cm-sized or 
larger dust aggregates to high number densities. In all 
such cases, knowledge about the low- and intermediate- 
velocity collision behavior of dust aggregates is of utmost 
importance to correctly describe the fate of the dusty 
components. In highly mass-loaded regions of PPDs, 
such as the dust sub-disk or instability regions, the coef- 
ficient of restitution determines the reduction of the rel- 
ative velocities among the dust aggregates. A low value 
of the coefficient of restitution eases the occurrence of 
gravitational instabilities. It is, thus, important to know 
the threshold velocity for the onset of fragmentation, so 
that the size evolution of the dust aggregates and, thus, 
the efficiency of the gravitat ional instability i n PPDs can 
be co rrectly determined (see lJohansen et all (|2006l 120081 

5.2. Debris Disks 

As relative velocities in virtually gas-free debris disk 
are typically larger than tens of meters per second, frag- 
mentation dominates the outcomes in collisions between 
dust aggregates. Thus, one would expect a broad size dis- 
tribution of dust aggregates from the monomer grains or 
the radiation-pressure blowout size (whatever is larger) 
to the largest occurring bodies. However, the recent 
discovery of ultra-cold debris disks with dust tempera- 
tures below the black-body equi librium temperatu re has 
severely challenged this picture (jEiroa et al.ll201lH . Such 
low temperatures require dust materials with very low 
absorption in the far infrared and the absence of a source 
for particles smaller than the wavelength. The former 
can possibly be reached with icy constituents, the lat- 
ter requires collision velocities below the fragmentation 
threshold. 

We found in our investigation that the fragmentation 
limit is as low as ~ 50 cm s _1 for large dust aggregates 
consisting of micrometer-sized Si02 grains. This is very 
close to the ~ 10 cm s" 1 sticking threshold for the 
monomer particles (jPoppe et alJl2000| ) . Unfortunately, it 
has only recently been p ossible to produce mic rometer- 
sized water-ice particles (|Gundlach et al.l 12011!) so that 
impact experiments with monomer particles and aggre- 
gates hav e not yet been p e rform ed. However, the experi- 
ments by lGundlach et all ()2011[ ) showed that the surface 
force for the water-ice particles is about a factor of ten 
higher than for silica so that one can expect a similar in- 
crease for the threshold velocities for monomer sticking 
and aggregate fragmentation, i.e. Uf ra g,ico ~ 10 m s _1 . 
Future experiments will have to show whether this is true 
or not. 

5.3. Planetary Rings 

In Saturn's main rings, the particle sizes have been es- 
timat ed to range from ~ 1 cm to ~ 10 m (|Zebker et al.1 
fl985h and are believed to be at least covered by a regolith 
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layer ([Poulet fc Cu zzi 2002), which makes them similar 
in their collision behavior to our dust aggregates. In the 
rings, the orbital shear leads to random velocities among 
the ring particles. This effect is counterbalanced by the 
energy loss due to the inelastic collisions among the ring 
particles, leading to a steady-state velocity distribution 
if and only if the coefficient of restitution decreases with 
increasing collision velocity. A basic description of this 
behavior is given in IGoldreich fc Tremaind lfl978h . The 
authors also found a direct correlation between the op- 
tical depth of the ring and the coefficient of restitution 
if the rings are in dynamical equilibrium. For reason- 
able optical depths for Saturn's main rings, the required 

coefficients of restitution are e ~ 0.6. Both conditions, 
the negative correlation between the coefficient of resti- 
tution and the collision velocity as well as the required 
values of the coefficient of restitution show that the phys- 
ical composition of the ring particles must be such that 
the threshold velocity between elastic and plastic im- 
pacts is just slightly lower than the typical impact speeds 

(~ 1 cm s _1 ). Clearly, much more work has to be done 
before we can really understand the collision behavior 
of ring particles (e.g. use icy aggregates, much lower 
collision velocities, include rotation, etc.), but the dust- 
aggregate collision experiments presented in this article 
show how to proceed. 

6. SUMMARY 

We investigated the low-velocity collision behavior of 
5-cm sized cylindrical dust agglomerates made by com- 
pression of micrometer-sized Si02 particles. We mea- 
sured the coefficient of restitution and the onset of frag- 
mentation for agglomerates with volume filling factors of 
<fi = 0.3 and <j> = 0.4 and found that (1) at the very lowest 
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velocities (v e g < 4 cm s _1 ), th e coefficient of rest i tution 
follows the solid state model bv lThornton fc Ninel (|1998f ) 
within the measurement uncertainties as shown by the er- 
ror bars in Figure|4] that (2) the coefficient of restitution 
decreases steeply to e = 0.2 for a volume filling factor of 
0.3 for velocities of 4 cm s _1 < v c fj < 6 cm s _1 and to 
e = 0.3 for a volume filling factor of 0.4 for velocities of 
4 cm s" 1 < v c ff < 10 cm s _1 ), that (3) the coefficient of 
restitution remains constant for higher velocities, until 
(4) fragmentation dominates for v c g > 40 ± 10 cm s _1 . 

Our own numerical simulation s, using the mo l ecular - 
dynamics approach described in iSeizinger et ail {2012), 
yield a reasonable match to the experimental data over 
the entire bouncing regime. We discussed the conse- 
quences of our results concerning PPDs, cold debris 
disks, and planetary rings. 
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